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Abstract 



It is found that the Hamiltonian of S=l/2 isotropic Heisenberg chain with sites and 
elliptic non-nearest-neighbor exchange is diagonalized in each sector of the Hilbert space 
with magnetization N/2 — M, 1 < M < [A^/2], by means of double quasiperiodic mero- 
morphic solutions to the M-particle quantum Calogero-Moser problem on a line. The 
spectrum and highest-weight states are determined by the solutions of the systems of 
transcendental equations of the Bethe-ansatz type which arise as restrictions to particle 
pseudomomenta. 
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In recent years, much attention has been paid to studies of ID lattice systems, due 
to their relevance to principal notions of field theory and experimental investigations of 
effectively low-dimensional crystals. Even the simplest lattice systems, namely isotropic 
S=l/2 Heisenberg chains, have unveiled rich structure and provided nontrivial examples 
of many-body interactions. The corresponding mathematical problem consists in finding 
the proper analytic tool for the diagonalization of the model Hamiltonian 

7YW = ^ E hU - k){a,a, - 1) h{j) = h{j + N) (1) 

i<j¥'k<N 

where aj are Pauli matrices acting on spin at jth site. 

At finite N, it has been successfully treated in the integrable cases of nearest-neighbor 
coupling solved by Bethe [1] 

h{j) = 5|i(modJV)|,l (2) 

and long-range trigonometric exchange proposed independently by Haldane and Shastry 
[2] 

M.j^f^^ta^^V'. (3) 



TT N. 

At present, a number of impressive results are known for both these models. In particular, 
they include the additivity of the spectrum under proper choice of "rapidity" variables 
[1,3], the description of underlying symmetry [4,5], construction of thermodynamics in the 
limit N ^ oo [6,3], the connection to the continuum integrable many-body problems [7,2], 
and closed-form expressions of correlations in the antiferromagnetic ground state. The 
rich collection of various generalizations and physical applications of Bethe and Haldane- 
Shastry models can be found in recent review papers [8,9]. 

Several years ago, I have introduced a more general one-parametric form of spin ex- 
change which provides another example of integrable lattice Hamiltonian (1) [10]. It 
has been motivated by the similarity of the Lax representation of the Heisenberg equa- 
tions of motion for continuum and lattice models. In the former case, the most general 
translationally-invariant integrable Hamiltonian with elliptic pairwise particle interaction 
has been found by Calogero [11] and Moser [12], 



HcM = 2 



e^+a(a+i)i:p(x, 

13=1 "-^P /3^7 



(4) 



The existence of extra integrals of motion commuting with (4) has been demonstrated in 
[13]. Recently, the eigenvalue problem for the elliptic Calogero-Moser operator received 
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much attention due to its relation to the representations of double affine algebras and 
solutions of Knizhnik- Zamolodchikov-Bernard equations [14,15]. 
The lattice analog of (4) is given by (1) with 





7r\2 r 


(- 


sm — 


Vvr 





PivO + -Civ - 



(5) 



where p^{x), Cjv('^) arc the Wcierstrass functions defined on the torus Tjv = C/ZiV + Za;, 
cu — in, K & R+. Remarkably, it turned out that the exchange (5) comprises both (2) and 
(3) [10]: in fact, the factor in (5) is chosen as to reproduce the nearest-neighbor couphng 
under periodic boundary conditions (2) in the limit k — > and the long-range exchange 
(3) in the limit k — > oo. 

However, till now much less is known about the lattice model with the exchange (5) 
in comparison with its limiting forms due to the mathematical complexities caused by 
the presence of the Weierstrass functions. The family of the operators which commute 
with Ti,^^^ has been found only recently [16]. The simpler case of infinite chain N oo, 
h{j) [sinh(7r//t)/ sinh(7rj'//t)]2 has been considered in detail in [17]. As for finite A^, the 
description of the spectrum has been performed only for simplest two- and three-magnon 
excitations over ferromagnetic vacuum [10, 18, 19]. 

The aim of this Letter is to demonstrate the remarkable correspondence between 
the highest-weight eigenstates of the lattice Hamiltonian with the elliptic exchange (5) 
and double quasiperiodic meromorphic eigenfunctions of the Calogero-Moser operator (4) 
which allows to formulate the equations of the Bethe-ansatz type for calculating the whole 
spectrum. 

The Hamiltonian (1) commutes with the operator of total spin S = | J2f=i (^j- Then 
the eigenproblem for it is decomposed into the problems in the subspaces formed by the 
common eigenvectors of S3 and such that S ^ S^^ N/2 - M, < M < [N/2], 

The eigenvectors > are written in the usual form 

TV M 



|V^(^)>= ^M(ni..nM) n ^nJO >, (7) 



where |0 >= | tt ••• T> is the ferromagnetic ground state with all spins up and the 

summation is taken over all combinations of integers {n} < N such that Vl/l<J<J^ii~'^'v) 7^ 
0. The substitution of (7) into (6) results in the lattice Schrodinger equation for completely 
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symmetric wave function ipM 



N M 

E E P^i'^P - s)'0M(rii, s, np+i, ..um) 

s^ni,..nM /3=1 



+ 



M 



i)M{ni, ..Um) = 0. 



(8) 



The eigenvalues {Em} are given by 



E 



M 



IT UJ 



J - sin - {Sm + - 



UJ 



2M{2M -1) - N 



UJ 



Civ - - MCi - 



UJ 



(9) 



where Ci(^) is the Weierstrass zeta function defined on the torus Ti = C/Z + Zuj. 
To find the solutions to (8), let us consider the following ansatz for ipM- 



(fiMini, ..Um) 



/ M \ 

cxp i-i ^PyuA Xm (ni, .mm), 



(10) 



where ttm is the group of all permutations {P} of the numbers from 1 to A^" and is 
the solution to the continuum quantum many-particle problem 



M 



0. 



;i2) 



It is specified up to a normalization factor by the particle pseudomomenta (pi, .-Pm)- The 
standard argumentation of the Floquet-Bloch theory shows that due to perodicity of the 
potential term in (3) Xm obeys the quasiperiodicity conditions [18] 



Xm{^i^ ■■^H + ..Xm) = exp{ipf3N)xM {xi, ..xm), 



(p)i 



(13) 



Xm (xi, ..x/3 + UJ, ..Xm) = exp(g^(p) + ip/3U})xMixi, -Xm), < ^m{qp) < 27r (14) 

1<P<M. 

The eigenvalue Em(p) is some symmetric function of (pi, ..Pm)- The set {q^ip)} is also 
completely determined by {p}. In this Letter I do not refer to the exphcit form of these 
functions which is still unknown for M > 3. 

The structure of the singularity of p]\[{x) Sit x — implies that Xm ^e presented 
in the form 

(15) 



Xm 



F^\x^,..Xm) ^, ^ n r ^ 

— ^, G[Xx, ..Xm) = 11 CTNiXa - Xp), 

U[Xi,..Xm) a^i3 
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where aN{x) is the Weierstrass sigma function on the torus T^. The only simple zero 
of crjv(a;) on is located Bit x — 0. Thus [G{xi, ..xm)]~^ absorbs all the singularities of 
Xm on the hypersurfaces Xa — xp. The numerator in (15) is analytic on (Tat)^ and 
obeys the equation 



M 



M 



2Em{p) - {pN{Xa - Xp) - C^i^a " X/^)) 



pip) 



(16) 



The regularity of the left-hand side of (16) as x^ — > x,y implies that 

f d d 



dx^ dx^^ 



FP{xi, ..xm] 







(17) 



for any pair (/x, v). 

The remarkable fact is that the properties (13-15,17) of Xm allow one to validate the 
ansatz (10-11) for the eigenfunctions of the lattice Schrodinger equation (8). Substitution 
of (10) to (8) yields 



M 



X! S X)'5/3(npi, ..npAf) + 

PGTTM 1/9=1 



M 



M 



A\npi,-npM)} =0, (18) 



where 



AT 



Si3(npi, ..upm) = pN{npp- s)Q^p^(pi\npi,..npM)- (19) 



s^npi,..npM 



The operator Q^^^ in (19) replaces /9th argument of the function of M variables to s. 

To calculate the sum (19), let us introduce, following the consideration of the hyper- 
bolic exchange in [17], the function of one complex variable x. 



M 



W'j^\x) = Y PN{npf3 -s- x)Q^f^^''^ (p%\npi, ..upm)- 

s=l 

As a consequence of (11), (13-14) it obeys the relations 

W\^\x + 1) = W^P{x), W^p^\x + uj) = e^^{qp{p))W'j^\x) 



(20) 



(21) 



The only singularity of Wj^^ on the torus Ti = C/Z -|- Zo; is located at the point x — Q. 
It arises from the terms in (20) with s — npi, ..npM- The Laurent decomposition of (20) 
near x — has the form 



W^^\x) = W^2X~^ + W^iX~^ + Wo + 0{x). 



(22) 
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The explicit expressions for w^i can be found from (20), 

W-2 = (fMinpi, ..npu) (23a) 

onpp 

+(-l)^G(ni, ..um) Tpxinpi, ..npM)Qf''^^ exp -i ^Pufip^ F^\npi, ..upm) 

X^0 \ v=\ ) 

(236) 
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Wq = Sfiinpx, ..Upm) + — (^^(npi, ..upm) + (-l)^G'(ni, ..hm) 
X J2 Tpxinpi, ..upm) [Upxinpi, ..npM)Qf'''^ + Piv(np^ - npx)dQf'''^] (23c) 

/ M 

X 



exp ^-i J^PuHp^^ F^\npi, ..upm), 



where 

rr r \ / \ TJ '^N{npp - npji) 

J-I3\\npi, ..npM) = cFNyripx - npp) [[ — -, 

p^/3,A (^N[npp - npx) 

Upx{npi, ..Upm) ^ pNi^PX - npp) - pN{npp - npx) CN{npp-npx), 

(— 1)^ means the parity of the permutation P and the action of the operator dQ^^^^^ on 
the function y of M variables is defined as 



dQ^;'''^Y{z,, ..Zm) = ^^(^1' ■■^M)U,=n,,. (24) 



The next step consists in writing the explicit expression for the function W^\x) obeying 
the relations (21) and (22) [17], 

Wl,^\x) = exp{apx) ^^l^^ ^ {w-2(pi(x) - pi{rpi) + {w-2{ap + 2Ci(r/3)) - W-i) 
ai[ri3 - X) 

x[Ci(x - rp) - Ci(x) + Ci(r^) - Ci(2r;3)]}. (25) 

The Weierstrass functions pi,Ci and ai in (25) are defined on the torus Ti and the 
parameters a^,r^ are chosen as to satisfy the conditions (21), 

= (7ri)-ig^(p)Ci(l/2) = -(47ri)-^g^(p). 

By expanding (25) in powers of x one can find Wq in terms of w_2,w_i,qf^ and obtain 
the explicit expression for iS^(npi, ..Upm) with the use of (23a-c). It turns out that the 
equation (18) can be recast in the form 



(p''P\npi, ..Upm) 
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= ^G{ni, .Mm) XI [^pxinpi, ■■npM) + Zxpinpi, ..upm)] , (26) 

where 

ff,(p) = (7ri)-^g^(p)Ci(l/2) - Ci((27ri)-^gMp)), (27) 
^/3(P) = ^Pi((27r0-^g^(p)) (28) 
and Zpx{npi, ..upm) is defined by the relation 

Zi3\{npi, .mpm) = Tfjx{npi, .mpm) Ufjx{npi, ..npM)Q'^^^^^ + pN{npx - npp) 

x(ag("^^) - fp{p)Q^;''^)\ exp {-lY^p.npA F^^\np,, .mpm). (29) 

Turning to the definition (11) of (/j'^^' one observes that each term of the left-hand side 
of (26) has the same structure as the left-hand side of the many-particle Schrodinger 
equation (12) and vanishes if £m and fpip) are chosen as 

fp{p) = -ipp, P=l,..M, (30) 

M 

Sm = Em(p) + J2 ^piP)- (31) 
1^=1 

Now let us prove that that the right-hand side of (26) also vanishes. The crucial observa- 
tion is that the sum over permutations in it can be recast in the form 

J2 (-1)^ H[^/3a(^pi> ■■npM) - ZxoinpRi, ..uprm)], 

where R is the transposition {j3 ^ A) which leaves other numbers from 1 to M unchanged. 
The term in square brackets is simplified drastically with the use of the identities 

TxisinpRi, --npRM) = Tpx{npi, ..urm), Ux/sinpRi, .mprm) = Upx{npi, ..upm) 

Qt^'^FinpR^, .mprm) = Q^;'"^F{np,, ..urm). 
Taking into account the relations (29-30), one finds 

Z,3x{npi, ..urm) - Zxp{npRi, ..urrm) = Tpxiripi, ..nRM)pN{npx - ripp) 



X exp 



M 

{Pp +Px)npx+ Pp'^Pp 



(32) 

The last factor in (32) vanishes due to the condition (17) imposed by the regularity of 
the left-hand side of the Schrodinger equation (16). 



It remains to show that the states of the spin lattice given by (7) with the functions ■^m 
of the form (10-11) are highest-weight states with S — S3. This statement is equivalent 
to the relation S+j^'^^^ >= 0, which can be rewritten as 

M 

EE E QfA\npi,..npM)^0, (33) 

13=1 p-=_(/3) S7^ni,..nM-i 

where {tt^m} are the subsets of ttm'- P G t^'m ^ — The sums over s in (33) can 
be reduced and presented in the closed form by using the technique described above. It 
turns out that the left-hand side of (33) contains the factors similar to the last factor in 
(32) and vanishes due to the condition (17). 

The descendant states with < S are obtained by acting with S_ on the basic states 
^ipW > (7). Thus the present consideration allows, in principle, to reproduce all the 
eigenvectors of Ti.^^^ for the exchange (5) as it has been done by Bethe [1] for nearest- 
neighbor spin coupling. The equations (30) for the pseudomomenta {p} constitute the 
analog of the usual Bethe ansatz. The spectrum is given by the relations (9) and (31). 

In conclusion, it is demonstrated that the procedure of the exact diagonalization of 
the lattice Hamiltonian with the non-nearest-neighbor elliptic exchange can be reduced 
in each sector of the Hilbert space with given magnetization to the construction of the 
special double quasiperiodic eigenfunctions of the many-particle Calogero-Moser problem 
on a continuous line. The equations of the Bethe-ansatz form appear very naturally as 
a set of restrictions to the particle pseudomomenta. The proof of this correspodence 
between lattice and continuum integrablc models is based only on analytic properties of 
the eigenfunctions. One can expect that the set of spin lattice states constructed by this 
way is complete. This is supported by exact analytic proof in the two-magnon case. 

The analysis of exphcit form of the equations (30) available for M — 2,3 shows that 
the spectrum of the lattice Hamiltonian with the exchange (5) is not additive being given 
in terms of pseudomomenta {p} or phases which parametrize the sets {p, q} [10,19]. The 
problem of finding appropriate set of parameters which gives the "separation" of the 
spectrum remains open. It would be also of interest to consider various limits {N 
00, K — > 0, 00) so as to recover the results of the papers [1,3,17] and prove the validity of 
the approximate methods of asymptotic Bethe ansatz after finding explicit form of the 
functions q/3{p) and Em{p)- 
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discussion. The support by the Ministry of Education, Science and Culture of Japan is 
gratefully acknowledged. 
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